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Abstract 

We present the results of a study of the azimuthal characteristics of 
ionospheric and seismic effects of the meteorite ’Chelyabinsk’, based on the 
data from the network of GPS receivers, coherent decameter radar EKB 
SuperDARN and network of seismic stations, located near the meteorite 
fall trajectory. It is shown, that 6-14 minutes after the bolide explosion, 
GPS network observed the cone-shaped wavefront of TIDs that is inter¬ 
preted as a ballistic acoustic wave. The typical TIDs propagation velocity 
were observed 661 ± 256m/s, which corresponds to the expected acoustic 
wave speed for 240km height. 

14 minutes after the bolide explosion, at distances of 200km we ob¬ 
served the emergence and propagation of a TID with spherical wavefront. 
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that is interpreted as gravitational mode of internal acoustic waves. The 
propagation velocity of this TID was 337 ± 89m/s which corresponds to 
the propagation velocity of these waves in similar situations. At EKB 
SuperDARN radar, we observed TIDs in the sector of azimuthal angles 
close to the perpendicular to the meteorite trajectory. The observed TID 
velocity (400 m/s) and azimuthal properties correlate well with the model 
of ballistic wave propagating at 120-140km altitude. 

It is shown, that the azimuthal distribution of the amplitude of ver¬ 
tical seismic oscillations can be described qualitatively by the model of 
vertical strike-slip rupture, propagating at Ikm/s along the meteorite fall 
trajectory to distance of about 40km. These parameters correspond to 
the direction and velocity of propagation of the ballistic wave peak by the 
ground. 

It is shown, that the model of ballistic wave caused by supersonic 
motion and burning of the meteorite in the upper atmosphere can sat¬ 
isfactorily explain the various azimuthal ionospheric effects, observed by 
the coherent decameter radar EKB SuperDARN, GPS-receivers network, 
as well as the azimuthal characteristics of seismic waves at large distances. 


1 Introduction 

The bolide, observed on February, 15 2013 above Chelyabinsk region, Russia, 
is the second powerful event since the fall of the Tunguska meteorite in 1908 
[Zotkin and Tsikulinl966llBen-Menaheml975| . According to [Popova et al.2013| , 
’Chelyabinsk’ meteorite had an initial size of about 19.8 meters and an ini¬ 
tial weight of about 13000 tons. When entering the atmosphere, its velocity 
was about 18.6 km/s. The most powerful explosion (the main explosive dis¬ 
ruption) occurred at 03:20:32UT (09:20:32LT) at « 27km height. After that, 
two more explosions occurred: second disruption (03:20:33.4UT) and third dis¬ 
ruption (03:20:34.TUT) [Popova et al.2013| .The energy of its burning in the 
atmosphere, that caused the destruction of buildings and windows, accord¬ 
ing to the calculations of [Popova et al.201^ was equivalent to 590 kilotons 
(kT) of TNT, while the part of the energy, associated with explosions, looks 
relatively small [Popova et al.2013 (see support materials(SM), pp.62-69 in 
[Popova et al.2013'| X The paper pilpatov et al.2013] has shown, that the ki¬ 
netic energy of the explosion of the ’Chelyabinsk’ meteorite was « 20kT based 
on the data of 11 infrasound stations of International Monitoring System of 
nuclear tests (MSM). 


Midscale travelling ionospheric disturbances (MSTlDs) with few hundred 
kilometers spatial scale, and few tens of minutes temporal scale, that were caused 
by the flight and explosion of the ’Chelyabinsk’ meteorite, were detected by vari- 


ous instruments [Alpatov et al. 20131 IGivishvili et al. 20131 

Gokhberg et al.2013 

[Berngardt et al.2013a| Yang et al.2014[ IR.uzhin et al.2014 

Chernogor2015 . 


By using low-orbital tomography at vertical ionospheric profiles over the 
European territory of Russia, wave-like disturbances (wavelength « 1.2° at lati¬ 
tude) were detected, that appear 2.5-3 hours after the explosion [Alpatov et al.2013 
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IGivishvili et al.2013) . The disturbances of the total electron content (TEC), 


caused by the explosion, were detected by ARTU GPS station Gokhberg et al.2013 
IRuzhin et al. 20141 |Chernogor2015| . According to [Gokhberg et al.2013| TEC 
variations had an inverted N-wave shape. The distribution of TEC variations 
amplitudes was not spherically symmetric. The ionospheric response to the ex¬ 
plosion had a number of features (waveform, axial symmetry in the ionospheric 
disturbances distribution) distinguishing it from the ionospheric response to 
earthquakes, tsunamis, explosions and volcanic eruptions [Gokhberg et al.2013 


In [Yang et al.2014[ , using GPS networks near the meteorite trajectory, in 
Japan and in the United States, three different types of MSTIDs were found: 
the higher-frequency (4.0-7.8mHz, periods 4.2-2.Imin) disturbances were ob¬ 
served around ARTU station and having « 862to/s velocity; the lower frequency 
(1.0-2.5mHz, periods 16.7-6.7min) disturbances, propagating with a mean ve¬ 
locity « 362to/s, were observed at distances of 300-1500km from Chelyabinsk; 
the higher-frequency (2.7-llmHz, periods 1.5-6.2min) disturbances with a mean 
propagation velocity « 733m/s were registered in USA region. Within 14min 
after the bolide explosion, on several nearby GPS-stations different TEC vari¬ 
ations were observed: some of them were caused by solar terminator and some 
of them propagate radially from the explosion epicenter to the distances up to 
500-600km with 320-360m/s velocity [Berngardt et al.2013a[ . At 06:00-10:00UT 
at most GPS stations, located away from the meteorite trajectory, an intensive 
TEC variations were observed, which had the form of wave packets with 30- 
40min period [Berngardt et al.2013a 


According to the over-the-horizon coherent EKB SuperDARN radar, located 
200 km from the meteorite explosion place, it is shown that the flight and fall 
of the meteorite caused the formation of MSTIDs with I50-200km horizontal 
spatial scale, propagating throughout the entire ionosphere (in E- and F-layers) 
to the distances up to 1500km from the epicenter. They have two typical prop¬ 
agation velocities - 220m/s and 400m/s [Berngardt et al.2013b[ . 

This paper presents the results of studying the azimuthal features of MSTIDs 
and seismic effects of ” Chelyabinsk” meteorite according to GPS network, EKB 
SuperDARN radar and network of seismic stations, located near the meteorite 
trajectory. 


2 Observations 

2.1 GPS TEC variations 

The study of TEC variations in the ionosphere during Chelyabinsk meteoroid 
flight was based on phase measurements of double-frequency GPS receivers 
(Figg. The data (with 1-sec sampling rate) of ARTU GPS-station included 
in the International GNSS Service (ICS) network were downloaded from the 
Scripps Orbit and Permanent Array Center website (http://sopac.ucsd.edu). 
GPS data (with 30-sec sampling rate) of stations TRIM and ORNB operated 
by NAVGEOKOM company were available from http://www.navgeocom.ru. 
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Data of GPS receivers in the Chelyabinsk region (CHEL, CH02, CH03, CH08, 
CH33, CH34, CH36, CH39; 5-sec sampling rate) were kindly provided by LLC 
’’GEOSalyut” (Moscow) and LLC’’Poleos” (Chelyabinsk). The Kazan GPS net¬ 
work data (stations LME, ATNI, CHEL, CHIS, KAYB, KZN2, LAIS, MAMA, 
MENZ, MKTM, NFTK, NOVO, OKTB, OPEN, SAMR, SARA, TYUL, UFAB, 
UTCH, ZAIN, ZDOL, 30-sec sampling rate; stations 3192, 3196, 3197, 3199, 
3204, 15-sec sampling rate) were provided by Federal Kazan University (KFU) 
and RPC ’’Geopoligon” KFU (Kazan). The last five GPS-stations are very close 
to each other and shown as single 3197 station at Fig{T] 

The area around F2-layer maximum (height hmax) make the main contribu¬ 
tion to TEC variations. Therefore, it is assumed that TEC variations are formed 
in an ionospheric piercing point (IPP), where IPP is the point of intersection 
of the ’’receiver-satellite” line of sight (LOS) with an ionospheric layer at hmax- 
According to ARTI ionosonde, located in Ekaterinburg, during meteorite fall 
hmax was about 243 km. We used this value in our calculations. 

To separate disturbances triggered by the meteoroid explosion, the initial 
TEC series I(t) were converted to an equivalent ’’vertical” value |Klobucharl986| . 
and then we filtered it in the period range 1-20 min [Afraimovich et al.l998] . 
TEC variations occurring on the day of the meteorite fall were compared with 
the TEC behavior observed on previous and subsequent days. The analyzed pe¬ 
riod featured quiet geomagnetic conditions: Kp-index did not exceed 1 within 
00:00-12:00 UT, and there were no strong solar flares. The seismic background 
was also quiet in the Chelyabinsk region and nearby areas. Some difficulties 
in separating the meteorite-generated ionospheric disturbances were caused by 
the fact that the meteoric impact occurred at sunrise, when the ionosphere 
exhibited strong variability. [Afraimovich et al.2009] demonstrated that solar 
terminator (ST) generates TEC wave disturbances in the ionosphere. Such dis¬ 
turbances can be seen for 2-4 h and even precede ST, which is attributed to the 
ST passage through the magnetoconjugate region. Between 02:00 and 06:00 UT 
on February 14-16, 2013, almost in all LOSs in the TEC variations there were 
quite intense oscillations with an amplitude of 0.2 TECU and a period of 15-17 
min. It is most likely that they were caused by the ST passage. 

Nevertheless, in spite of the presence of ST-caused TEC oscillations, in some 
receiver-satellite LOSs we managed to separate disturbances with a characteris¬ 
tic shape corresponding to the shape of a shock acoustic wave. In Figj^-h, such 
disturbances /STEC are marked with dashed lines. These disturbances were os¬ 
cillations with a period of « 10 min and amplitude 0.07-0.5 TECU exceeding 
the level of background oscillations on the reference days. 

To draw a spatial-temporal picture of the TEC response to the meteoroid 
explosion, we have constructed TEC disturbance maps (Figj^-f). Figure 2a- 
f presents spatial distribution of minima (blue dots) and maxima (red dots) 
of TEC oscillations (at hmax = 243A;to) over a period 03:26-04:00 UT. Thick 
concentric lines indicate the position of the disturbance wave front which can 
be traced through TEC variations minima (blue) and maxima (red). The first 
TEC disturbances were registered 6 min after the meteorite explosion at dis¬ 
tances of 80-100 km from the airburst location. Over a period 03:26-03:33 UT, 
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the TEC disturbances are characterized by fast dynamics. The front of distur¬ 
bances has a form of the cone dispersing from the meteorite trajectory. The 
angle at a cone peak changes from 4.9° to 22.9° (the average angle value is 
about 14.3°). The average velocity was estimated from movement of front lines: 
V = 661±256m/s. Since 03:34 UT the behavior of the TEC disturbances corre¬ 
sponds to the spherical wave front propagating from the center with coordinates 
54.90°A^; 60.85°if (blue cross in Figj^-f). The wave center is shifted approx¬ 
imately 37 km northwestward of the airburst. It coincides with the location 
of meteorite tertiary disruption (54.905°iV; 60.865°if) [Popova et al.201^ . The 
analysis suggests that the TEC disturbances propagated radially up to 600-700 
km and had wavelength around 220 km. The average velocity calculated from 
movement of front lines was V = 337±89m/s. It should be noted that spherical 
wave is generally registered eastward, northeastward, and southeastward from 
the airburst. The interference of two wave types (spherical and cone-shaped) is 
observed westward, northwestward, and southwestward from the airburst. This 
effect is more expressed during 03:34-03:38 UT. 


2.2 Coherent radar observations 


Measurements by GPS network are known to have difficulty in separating the 
spatial and temporal variation of ionospheric parameters. Unlike GPS data, 

SuperDARN radars [Chisham et al2007] have high spatial and temporal reso¬ 
lution, and allows us to differ spatial and temporal effects. The first Russian 
coherent radar EKB is located 56.5°fV, 58.5°E (near Ekaterinburg city), approx¬ 
imately 200 km to the North-West from the meteorite fall region and operated 

with 60km spatial and Imin temporal resolution in the sector —6°-44° from 

North at distances up to 3500km. Fig{^ shows the location of the radar and the 
orientation of each of the 16 rays (beams) of its field of view. 

After the meteorite fall, EKB radar recorded a large variety of ionospheric ef¬ 
fects [Berngardt et al.20I3a[|Berngardt et al.20I3bl|Berngardt20I3c]|Kutelev and Berngardt20~r^ . 
The most powerful and basic effect is the observations of formation and ra¬ 
dial passage of MSTIDs at E- and F-layer heights [Kutelev and Berngardt20I3| . 

Fig[^ presents the power of the scattered signal according to EKB radar data at 
3,6,9,12 beams (Figj^-L) during February 15, 2013 02:00-05:00UT. In addition, 
the figure shows the power of the scattered signal in geomagnetically quiet days 
- February 09, II, 18, 2013. The choice of these days been investigated in detail, 
for example, in Kutelev and Berngardt20T3] . 

The tracks marked by areas 1-2 in the figure correspond to refraction of the 
radiosignal at ionospheric irregularities, which propagate radially from the radar 
with velocities 400 and 200 m/s at E- and F-layer altitudes [Kutelev and Berngardt20I3| . 

According to Figj^ the ionospheric effects in the field of view of the radar 
have strong azimuthal dependence. In particular, they present on 6 beam, and 
do not appear at 3 and 12 beams. To emphasize the azimuthal effects, we used 
the following technique: 

- For each moment, for each beam, and for each range (this allows us to 
uniquely define each region in the field of view with approximately lOOxIOOkm 
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size) we calculate the average power (in dB) higher than the noise level OdB, 
over the magnetically quiet days February, 09,11,18, 2013; 

- Then, when processing 02/15/2013 data for each moment, beam and range 
we analyse only the data exceeding the threshold level (tripled the calculated 
average power at this point). 

- The power and the number of observations, accumulated for the given time 
period [Tstart,Tstop] for every fixed beam and every fixed range, are shown in 
Figll) 

Fig[^,C,E shows the intensity of the scattered signals in the radar field 
of view above the threshold level. Fig[^,D,F presents the number of signals 
exceeding a threshold level. 

As one can see from Fig|^-D, the main effect is concentrated in the sector 
of directions 2° — 18° from the north. The dotted line marks the direction 
corresponding to the perpendicular to the meteorite trajectory. 

2.3 Seismic observations 

For the analysis of the seismic characteristics, we used seismograms obtained by 
worldwide networks of the broadband seismic stations Iris/Ida (II, http://ida.ucsd.edu/) 
and Iris/USGS (lU, http://earthquake.usgs.gov/regional/asl/) and by four re¬ 
gional networks: Iris/China (IC), Kazakhstan (KZ), Kyrgyzstan (KR) and 
Baikal regional seismological center (BY, Russia). Analysis of seismograms is a 
bit complicated due to presence of seismic waves from Tonga earthquake (Febru¬ 
ary 15, 2013, 03:02:23UT, magnitude 5.8, coordinates 19.72S, 174.48W, depth 
71.6 km). In the rest, the seismic situation was very quiet - no local or regional 
earthquakes were registered on the seismograms. Visual analysis revealed the 
presence of the surface seismic waves at some stations that may be associated 
with meteorite according to their arrival times. These waves were recorded at 32 
seismic stations (Figj^), located at distances tq from 252 to 3654 km from the 
epicenter, and they are short-period oscillations with a period T = 3-16 sec and 
with duration up to 1 minute. For stations situated more than 3600 km from the 
fall site, a surface wave is lost in the seismic waves from the Tonga earthquake. 

The group velocity Cr of Rayleigh waves for periods of 20-100 seconds varies 
from 2.3 to 3.5 km/s. 

For the analysis of azimuthal distribution of surface seismic waves the maxi¬ 
mal amplitudes of vertical oscillations (Figj^), and the spectral corner fre¬ 
quencies fc were analyzed (Figj^). The distribution of seismic stations around 
the point of the explosion is non-uniform: there is a very big azimuthal gap ap¬ 
proximately from 330° to 70°. For this azimuthal range it is impossible to 
allocate the surface wave from the records due to the superposition of waves 
from the Tonga earthquake, as well as to the lack of broadband seismic stations 
( Fig[^). Analysis of at the nearest stations (tq < 620 km) shows that 
the maximum = 1414.4 nm is observed at station BRVK (ro= 616 km, 
azimuth - 104°), at the same time, the amplitude at station closer to the fall 
point is smaller: at ARU station (rp = 252 km, azimuth 317°) =1177.5 nm, 

at ABKAR station ( ro = 618 km, azimuth 189°) - = 762.4 nm (Figj^). 
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The azimuthal distribution of the maximum amplitudes of the Rayleigh 
waves at different stations indicates the absence of spherical symmetry (Figj^) 
- the maximum amplitudes are oriented along the direction of the meteorite 
trajectory, the minimum amplitudes are observed in perpendicular directions. 
There is three exceptions: an abnormally low amplitudes are observed at LSA 
(distance tq = 3654 km) and MKAR (rg = 1709 km) stations, located in az¬ 
imuths 113° — 127°, and abnormally high amplitude at the KONO station (az¬ 
imuth 301°, To = 3073 km). But statistically, the azimuthal distribution of the 
Rayleigh waves amplitudes has two peaks elongated with the meteorite tra¬ 
jectory: the first maximum is oriented towards the meteorite fall direction and 
the second - in the opposite direction. 

The analysis of the azimuthal distribution of showed that the distribution 
elongated with the meteorite fall trajectory (Az = 97°). At the same time, we 
found no spherical symmetry in seismic wave radiation that corresponds to 
bolide explosions. This allows to suggest that the seismic effect is caused by 
single extended event rather than a series of explosions. This fits well in with 
the model [Popova et al.2013] . An explanation of the observed asymmetry in 
the seismic wave amplitudes from the single point source is also not satisfactory. 
In some cases, azimuthal nonuniformity of the seismic wave amplitudes may be 
explained by an anisotropy of medium. But this explanation does not work 
in our case, as the main tectonic structures in the Ural region are oriented 
orthogonally to the meteorite trajectory. Thus, we can conclude that there is a 
directivity of seismic radiation. 

The effect of directivity is reflected both in the seismic wave amplitudes and 
on the corner frequencies. The concept of the corner frequency is associated 
with the rupture parameters (its velocity and length) and directivity function 
|Ben-Menaheml961| : 


Cr/V + cose sin [^^{Cr/V -cose)) 
CrlV-cosO sin {Cr/V + 0030 )) 


( 1 ) 


where Cr is the seismic wave velocity, V is the rupture speed, b is the length of 
rupture, uj is angular frequency, 9 is the angle between the rupture propagation 
direction and the direction to the seismic station. 

This function has a series of maxima and minima at frequencies, their posi¬ 
tions depend on the values h and V. The first maximum/minimum corresponds 
to a frequency fc = w/(27r) called the corner frequency. Thus, the corner 
frequency fc includes the effect of the propagation of rupture and it may be 
measured directly from the surface wave spectrum. It corresponds to the fre¬ 
quency of the inflection point of the seismic spectrum, separating the horizontal 
plateau at low frequencies and decaying oscillations with sloping envelope at 
high frequencies. For the meteorite surface waves, the corner frequencies vary 
widely (0.06-0.43 Hz), but the maximum number of stations observe frequencies 
0.15 - 0.25 Hz (Fig(^). 
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3 Models 

3.1 Acoustic signal model 

To explain the azimuthal characteristics of ionospheric effects, we made math¬ 
ematical simulations of acoustic signal propagation at different heights. 

Models based on one or more explosions are often considered as model 
caused by the ionospheric disturbances. However, from our point of view, the 
model of ballistic wave described in (SM in [Popova et al.201^ , pp.62-69) is 
the most realistic one. It successfully explains the observed destructions (which 
are associated with a passage of powerful acoustic wave) as continuous for¬ 
mation of perturbations along the meteorite trajectory. Similar models have 
been successfully used for explaining the zone of destruction after the Tun- 
guska meteorite fall [Zotkin and Tsikulinl966l IBen-Menaheml975] , which indi¬ 
cate a high degree of this model reliability. In this model, the amplitude of 
the source along the trajectory is determined by the luminosity of the bolide 
[Zotkin and Tsikulinl966| . The disadvantages of using other models to describe 
the destruction zone in the case of ’Chelyabinsk’ meteorite were investigated in 
detail in (SM in [Popova et al.20I3| ,pp. 62-69). 

Therefore, we used the simple model of the summation of partial pressures 
coming from each point of the meteorite trajectory with taking into account 
the acoustic propagation delay. We used the model of a moving source. The 
initial acoustic pressure disturbance is proportional to the luminosity of the 
bolide at the meteorite position for given moment [Zotkin and TsikulinI966l 
IBen-MenahemigTsl . The luminosity is shown in (SM in [Popova et al.20I3|, 
pp26-27). 

As a result of interference, the meteorite fight forms the focus of the acoustic 
wave - Mach cone (ballistic wave). In a homogeneous medium for the meteorite 
velocity 17.5km/s this wave spreads nearly perpendicularly to the flight trajec¬ 
tory. In a medium with nonuniform acoustic speed height profile the Mach cone 
shape substantially transforms. This transformation is determined by many 
factors, including the acoustic speed at different heights. Therefore at different 
distances, the angle between propagation direction and the meteorite trajectory 
becomes different. 

For estimates, we used the qualitative model for dependence of the acoustic 
wave pressure on the neutral density and acoustic speed [Yuenl96^ : 

V p{^o)Cs{Ro) 

Here, the delay Tg is determined by trajectory of the acoustic signal and 
based on the geometrical optics laws for the acoustical signal in inhomogeneous 
media [BlokhintsevI952[ [Clay and MedwinI977[ , where the refractive index is 
determined by the acoustic speed C's(^). The p( A) is the density of the neutral 
atmosphere and Ps{RQ,t) = aI{Ro,t) is the initial acoustic pressure of the 
source, which (in the first approximation) is equal to the luminosity of the 
























bolide I{Ro, t) at this point at corresponding moment |Zotkin and Tsikulinl966l 
IBen-Menaheml975] with accuracy of a constant multiplier a. Rq = Ro{t) is the 
position of the bolide in the moment t. 

The resulting partial acoustic pressure Pg at each point of {^,t) is summed 
up with taking into account the geometroacoustic propagation delay Tg from the 
bolide (source) point to the specified point. Since the source of the acoustic wave 
at each moment is assumed isotropic, we conducted an additional integration 
over all the acoustic signal propagation directions from the source point. 

The duration of the meteorite fall in the upper atmosphere was « 16sec, 
which corresponds to a wavelength « 5km, considerably less than a hundred 
kilometers - the distance to the studied effects. This can qualitatively substan¬ 
tiate the applicability of geometrical acoustics for qualitative estimates of the 
observed ionospheric effects. In the calculations of the acoustic field, we as¬ 
sumed that the acoustic speed and the neutral density profile are functions of 
the height H. The earth surface is considered to be flat, due to the smallness of 
the investigated distances. The acoustic speed and neutral density profile were 
estimated from the reference model NRLMSIS-00 for observation point. 

As an initial pressure perturbation Pg{Ro,t), forming an acoustic signal we 
used a short pulse propagating at the velocity Vq = 17.5fcm/s along the mete¬ 
orite trajectory from 90km to 25km altitude. Pg{Ro,t) depends on the height 
and luminosity curve defined in (SM in [Popova et al. 20131 , pp26-27). We used 
the data given in [Borovicka et al.2018] for entry point tq, the trajectory and 
meteorite fall velocity Vo¬ 
lt should be noted that the model is very rough, and does not include nei¬ 
ther details of the acoustic signal amplitude changes with height, nor conversion 
of neutral density variations into electron density variations. The effects also 
depends on the height [Maruyama and Shinagawa2014| and their presence is re¬ 
quired for a correct estimation of the amplitude of ionospheric effects. However, 
as it will be shown below, the model is suitable for qualitative explanation of 
azimuthal and temporal effects. 

The evolution of amplitude at 120km altitude, in E-layer, as a function of 
time, is shown at Figj^-c. One can see from the figure that the main acoustic 
effects should be observed 4-6 minutes after the explosion and spread further 
at the acoustic speed (« 400m/s) in a direction perpendicular to the meteorite 
fall trajectory. 

Fig[§l -f presents the evolution of acoustic wave pressure at 240km, in the 
F-layer, as a function of time. It shows that the main acoustic effects should 
be observed 7-8 minutes after the explosion and should propagate with acoustic 
speed (« 800m/s) in a direction of 40° — 50° to the meteorite trajectory. 

Fig[7}r-c shows the modeled acoustic wave maximal pressure at various heights 
and at various distances from the main bolide explosion site, with taking into 
account the acoustic sound speed prohle at the upper atmosphere. 

It can be seen that the acoustic wave pressure at different altitudes has a 
strong azimuthal dependence. At relatively low altitudes (below 150 km), the 
highest acoustic wave pressure is predicted at azimuths perpendicular to the 
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meteorite trajectory. This fits well with the results of [Popova et al.201^ which 
associate the azimuthal features of the destruction zone on the earth’s surface 
with an acoustic wave pressure. This justifies the possibility of using the simple 
model (H) for estimation of azimuthal effects. At altitudes above 150 km, the 
model predicts an increase of the acoustic wave pressure at azimuths « 45° to 
the meteorite trajectory in the direction of the meteorite fall. 

From Figj^-C one can see that the highest acoustic wave pressure is ob¬ 
served under the trajectory of the meteorite. Figj^ shows the dynamics of the 
acoustic wave pressure along the meteorite trajectory projection to the Earth’s 
surface. As it is shown at Fig0F, three acoustic disturbances propagate along 
the surface of the Earth: the most powerful propagates back towards its fall 
with « 1km/s velocity, which is defined by the meteorite trajectory fall angle, 
the meteorite speed and the acoustic speed profile. This velocity roughly coin¬ 
cides with the estimates given in (SM in [Popova et al.201^ , p.70), which also 
confirms the possibility of using our model to estimate the effects. The two 
other waves are acoustic ones, having velocity « 320to/s, which is associated 
with the propagation of spherical acoustic wave at ground altitudes. 

3.2 Seismic signal model 

The propagation velocity of surface seismic waves was about 3-3.5km/s [Tauzin et al.201,^ . 
Therefore the seismic oscillations at longer distances can not be directly caused 
by atmospheric acoustic waves having the velocity 10 times smaller. However, 
the movement of a supersonic powerful acoustic pulse (ballistic wave) on the 
ground can cause the observed seismic effect indirectly. This model, for ex¬ 
ample, was used for an explanation of seismic effects of Tunguska meteorite 
|Ben-Menaheml975] . To estimate the characteristics of the effect we carried 
out numerical calculations of the signal intensity in a simple model developed 
for the analysis of seismic signals generated during formation of geological faults. 

As the fault line we use the trajectory of the acoustic signal on the ground. Pa¬ 
rameters of fault formation (length, velocity and direction) are equal to the 
parameters of the acoustic pulse movement on the ground. 

To calculate the theoretical seismic radiation pattern we used the model 
proposed in [Ben-Menaheml96T| for a linear fault. According to this approach, 
the azimuthal distribution of the seismic waves is defined by the type of the 
fault and by the parameters of its formation: the rupture speed P, its length 6, 
the frequency of seismic waves w, and the phase velocity of Rayleigh waves Cr- 

In our case, two types of ruptures with the vertical position of the plane 
displacement may be considered as a seismic analogue - vertical strike-slip fault 
and the vertical dip-slip fault. In both cases, the maximal amplitudes have to 
orient according to the slip vector: for the strike-slip fault they are oriented in 
the direction strike fault, and for dip-slip fault they are perpendicular to the 
fault plane |KasaharaI9811 lUdfas et al2014] . 

In our analysis, we used the following formulas that qualitatively describe 
(up to a constant multiplier) the azimuthal distribution of amplitudes of vertical 
seismic oscillations in cases of different fault types. Thus, the amplitude of 
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the seismic waves depends strongly on the angle 6 between fault propagation 
direction and the direction to the seismic station. The amplitude also depends 
on the distance to the observation point rg . 

According to [Ben-Menaheml96T| , in the case of vertical strike-slip fault the 
azimuthal distribution of the Rayleigh displacement for the vertical component 
is given by: 


where 


ut 


cosO sinXr 


Xr 


Lob 



(3) 

(4) 


In the case of vertical dip-slip fault the azimuthal distribution of the Rayleigh 
displacement for the vertical component is given by: 


where 


1 sinYr 


y, = 


ujbcosO 

2Cr 


(5) 

( 6 ) 


For the simulation we choose the propagation line of the fault as the trajec¬ 
tory of the meteorite fall (azimuth of propagation 97°). 

We suppose the rupture speed V « Iktrijs, its length b « 40fcm, the fre¬ 
quency of seismic waves w, defined from fig[^ and the phase velocity of Rayleigh 
waves Cr ~ 3Akm/s. 

The direction of the rupture was opposite to the direction of the fall. This 
corresponds to the direction and the trajectory of the ballistic wave peak on the 
ground. 


4 Discussion 

Fig[^,B presents a comparison between observed MSTIDs and acoustic wave 
pressure distribution at 120km. One can see that the model of acoustic waves in 
the lower ionosphere describes the experimental data sufficiently well. The same 
conclusion can be made from a comparison of Figj^ and Fig[^ According to 
the estimates carried out before, the velocity of azimuth-dependent ionospheric 
irregularities was about 400m/s [Kutelev and Berngardt2013| . This corresponds 
well to the velocity of acoustic wave pressure variations at altitudes 120-140km. 
So we can conclude that the MSTIDs, observed by EKB radar, are well described 
by the effect of the ballistic wave passage at altitudes of 120-140km. This also 
stimulates similar disturbances in the charged components of ionospheric plasma 
(for example, by mechanism described in [Maruyama and Shinagawa2014| ). 

As simulations show, the acoustic radiation pattern at 240km has a specific 
two-horned structure, which describes the experimental results observed in the 
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first 14 minutes after the explosion sufficiently well. Fig|^-D shows the impo¬ 
sition of model calculations and experimental observations. We can see from 
Fig[^ that not only irregularities azimuthal distribution, but also the expected 
moment of irregularities observation are in a good agreement with the model. 
Estimates of the propagation velocity of TEC wave cone 6-14min after explosion 
is 661to/s ± 256m/s, and are in a good agreement with the calculated acoustic 
wave pressure velocity at 240km height (Figj^E). The observation of such waves 
by the EKB radar is hampered by the fact that the azimuthal effect is located 
outside the antenna pattern. 

14 minutes after the explosion we observed a formation and propagation of 
a spherical MSTID with velocity 337 ± 89m/s. Similar midscale (wavelength of 
200-300 km) concentric MSTIDs with velocities 138 —423m/s were detected by 
GPS receivers after the Tohoku earthquake (March 11, 2011) [Tsugawa et al.2011| 
The earthquake-generated MSTIDs began to be observed at distances « 300 A:to 
from the epicenter, and propagated up to lOOOkm. According to the numerical 
simulation |Matsumura et al.2011| , these propagating spherical MSTIDs can be 
related with gravitational mode of internal atmospheric waves. 

Thus, by GPS network we observed the two MSTID kinds in TEG: the con¬ 
ical wavefronts associated with ballistic wave 6-14 minutes after the explosion, 
and the spherical wavefronts associated with internal airwave 14 minutes after 
the explosion. 

Fig[^ shows the azimuthal distribution of the amplitude of vertical seismic 
oscillations, multiplied by y/ro and the model azimuthal distributions ([^[^. The 
model azimuthal distribution was calculated as the average of azimuthal distri¬ 
butions of expected vertical seismic oscillations, with taking into account the 
experimental distribution of their measured frequency (see Figj^D). Azimuthal 
distribution of normalized amplitudes, observed experimentally, is shown with 
blue crosses, the expected model azimuthal distributions are shown by lines. 
Green line corresponds to the case of vertical dip-slip fault, black and red lines 
correspond to the case of the vertical strike-slip fault. The black and red corre¬ 
spond to the cases of fault propagation in opposite and direct directions of the 
meteorite fall respectively. 

The figure shows that the model of vertical strike-slip fault describes qual¬ 
itatively the observed azimuthal distribution of the amplitudes of the vertical 
oscillations, and its elongation with the trajectory of the meteorite. This allows 
us to interpret ballistic wave peak, moving by the ground with Ikm/s velocity, 
as an equivalent vertical strike-slip fault, generating seismic waves, propagating 
with Rayleigh velocity. 


5 Conclusion 

We present the results of a studying the azimuthal characteristics of ionospheric 
and seismic effects of the meteorite ’Ghelyabinsk’ based on the data from the 
GPS receivers network, coherent decameter radar EKB SuperDARN and net¬ 
work of seismic stations located near the meteorite fall trajectory. 
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It is shown that 6-14 minutes after the bolide explosion network of GPS 
receivers steadily observed the cone-shaped wavefront of MSTIDs that is inter¬ 
preted as a ballistic acoustic wave. We assume the 240km effective ionospheric 
height. The typical MSTIDs propagation velocity were observed 661 ± 256m/s, 
which corresponds to the expected acoustic wave speed for 240km height. 

14 minutes after the bolide explosion, at 200km distance we observed the 
emergence and propagation of a MSTID with spherical wavefront, that may be 
interpreted as gravitational mode of internal acoustic waves. The propagation 
velocity of this MSTID were 337±89to/s which corresponds to the propagation 
velocity of these waves in similar situations. 

At EKB SuperDARN radar, we observed MSTIDs in the sector of azimuthal 
angles close to the perpendicular to the meteorite trajectory. The observed 
MSTID velocity (400 m/s) and azimuthal properties correlate well with the 
model of ballistic wave propagating at 120-140km. 

Thus, our modeling shows that ballistic wave produced by supersonic motion 
of the meteorite, and its burning in the upper atmosphere can satisfactorily ex¬ 
plain the various azimuthal ionospheric effects occurring during the first minutes 
after the explosion. The azimuthal dependence of ionospheric effects is observed 
by the coherent EKB SuperDARN radar and GPS network. 14 minutes after 
the explosion in the ionosphere additional spherically symmetric perturbations 
are beginning to occur that are associated with generation of gravitational mode 
of internal atmospheric waves. 

It is shown that the azimuthal dependence of the amplitude of vertical seis¬ 
mic oscillations can be described qualitatively by the model of vertical strike-slip 
fault, propagating at the velocity Ikm/s along the trajectory of the meteorite 
to the distance of about 40km. These parameters correspond to the velocity of 
ballistic wave pulse propagation by the earth surface. This allows us to consider 
the ballistic wave to be a moving source of seismic oscillations for a qualitative 
description of seismic observations during the meteorite fall. 

Our modeling and experimental data show that ballistic wave, caused by the 
supersonic flight of the bolide, plays a decisive role for the observed azimuthal 
dependence of ionospheric and seismic effects. 
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Figure 1: Measurements scheme during the meteorite ’Chelyabinsk’ 

fall. The gray dashed line marks the position of the solar terminator for 
03:20:32UT. Thick straight line shows the meteorite trajectory according to 
[Popova et al.2013| , the cross marks the strongest explosion. 
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Figure 2: Propagation of the TEC disturbances (a-f). The meteorite trajectory 
is shown according to [Popova et al.20^]^ by brown straight line. Three main 
disruption episodes are marked by brown crosses as given in [Popova et al.20T3] . 
The largest brown cross displays the position of main explosive disruption (the 
airburst). Red dots corresponds to TEC maxima positions, blue dots corre¬ 
sponds to TEC minima positions. TEC disturbances registered at GPS-station 
SAMR (satellite PRN26) (g), CH39 (satellite PRN06) (h). The airburst time is 
marked by a vertical dashed line. 
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Figure 3: The intensity of the scattered signal during 02:00-05:00UT at 3,6,12 
beams during geomagnetically quiet days February 09,11,15 and 18, 2013. Zone 
at H) marks the effect of refraction on MSTlDs arising after the flight and the 
explosion of the meteorite. 
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Figure 4: The normalized average scattering cross-section (A, C, E) and the 
number of ionospheric signals (B, D, F) in the field of view of the radar February 
15, 2013, that exceed the threshold level for the appropriate geographic area, 
calculated over the geomagnetically quiet days. Estimates of MSTIDs are made 
for the periods: A-B) - 01:30-03:00 (before the meteorite fall); C-D) - 03:00- 
04:30 (during the meteorite fall); E-F) - 04:30-06:00 (after the meteorite fall); 
G-H) - total number of ionospheric signals (G), and the scattering cross section 
(H), as a function of sounding azimuth^^rom the north direction). The vertical 
dashed line is the perpendicular to the meteorite trajectory. 
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Figure 5: A) The azimuthal distribution of seismic stations (triangles) relative 
to the point of the main bolide explosion (the center). The line marks the trajec¬ 
tory of the meteorite fall. B) The azimuthal distribution of vertical amplitudes 
reduced to epicentral distance 100 km. Dots show all the observed data (except 
data from ARU, BRVK, ABKAR stations with epicentral distance i620 km). 
C) The azimuthal distribution of amplitudes of vertical seismic oscillations at 
stations ARU, BRVK, ABKAR (shown by triangles), gray line marks the tra¬ 
jectory of the meteorite fall, the center corresponds to the point of the main 
bolide explosion. D) The distribution of vertical oscillations frequency over the 
seismic stations. 
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Figure 6: The dynamics of the acoustic wave pressure at 120km (left) and 240 
km (right). Arrow marks the meteorite trajectory. 


23 







Delay (min.) 





300m/s 


1000 m/s 


0 2 4 6 8 10 12 14 

Delay (min.) 


120 

110 

100 

90 

80 

70 

60 

50 

40 


Figure 7: Model of acoustic wave pressure at different heights: A) - 240km; B) 
- 120km; C) - 0km; D) - azimuth (from the north) of ballistic wave propagation 
at 240 km and 120 km heights; E) - velocity of ballistic wave propagation at 
240 km and 120 km heights; F) - propa^tion velocity of ballistic wave maximal 
intensity on the Earth surface. 
















































































Figure 8: A-B): Comparison of scattered signal power variations that are ob¬ 

served by EKB radar (A) and the model estimates of the acoustic wave pressure 
at 120 km height (B); C-D): TEC disturbances and distribution of acoustic wave 
pressure. Blue dots correspond to minimums in AT EC, red ones correspond to 
maxima in AT EC. Lines mark approximations of the wavefront over the data. 
Additionally shown C) - a model distribution of the acoustic wave pressure at 
03:30UT at 240 km height, D) - a model of the acoustic wave pressure for the 
period of 03:20-03:50UT at 240 km. E) The azimuthal distribution of the am¬ 
plitudes of vertical seismic oscillationgg multiplied by yTo (blue crosses) and 
the model azimuthal distribution of amplitudes of vertical seismic oscillations 
(Hi (lines). The expected distribution in the case of vertical dip-slip fault is 
shown by green line, black and red lines show the azimuthal distribution of the 
amplitudes for the vertical strike-slip fault. The black and red lines correspond 
to the cases of opposite and direct fault propagation relatively to the meteorite 
fall direction. 










